Mon. Not. R. Aston. Soc. 000.1111171(2002') Printed 2 February 2008 (MN IATpX style file v2.2) 



X-ray Preionisation Powered by Accretion on the First Black Holes. 
I: a Model for the WMAP Polarisation Measurement 



o 
o 

(N 



Massimo Ricotti and Jeremiah P. Ostriker 

Institute of Astronomy, Madingley Road, Cambridge CB3 0HA 
ricotti @ ast. cam. ac. uk, jpo @ ast. cam. ac. uk 



Accepted — . Received — ; in original form 3 1 October 2003 



ABSTRACT 

In this paper we investigate the possibility that there is a first phase of partial ionisation due 
to X-rays produced by black hole accretion in small-mass galaxies at redshifts 7 < z < 
20. This is followed by complete reionisation by stellar sources at z ~ 7. This scenario is 
motivated by the large optical depth to Thompson scattering, r c ~ 0.17 ± 0.04, recently 
, measured by WMAP. But it is also consistent with the observed Gunn-Peterson trough in the 

J> ■ spectra of quasars at z ~ 5 — 6. We use a semianalytic code to explore models with different 

C*~) \ black hole accretion histories and cosmological parameters. We find that "preionisation" by 

. X-rays can increase the intergalactic medium (IGM) optical depth from t c k 0.06 given 

by stellar sources only, to 0.1 ^ t c H 0.2, if a fraction of baryons 10~ 5 ;$ lo &c £ 10~ 4 is 
accreted onto seed black holes produced in the collapse of low metallicity, high mass stars 
before z ~ 15. To be effective, preionisation requires a non-negligible star formation in the 
first small-mass galaxies in which seed black holes are formed. By z ~ 20 — 25 the IGM 
is re-heated to 10, 000 K and the ionisation fraction is about 20%. The increase of the IGM 
• Jeans mass is effective in reducing star formation in the smaller-mass haloes. Large values of 

t are obtained in models with top-heavy stellar initial mass function only if pair-instability 
supernovae are not important. Seed black holes are assumed to accrete at near the Eddington 
£^ ■ limit with a duty cycle that decreases slowly with increasing time. Alternatively, a moderate 

^ \ fraction of the black holes must be ejected from the host galaxy or exist without merging into 

& ■ the supermassive black holes in galactic centres. The model predicts that dwarf spheroidal 

galaxies, if they are preserved fossils of the first galaxies, may host a mass in black holes that 
is 5-40% of their stellar mass. The redshifted X-ray background produced by this early epoch 
■ of black hole accretion constitutes about 5 — 10% of the X-ray background in the 2-50 keV 

5— i ' bands and roughtly half of the currently estimated black hole mass density was formed at early 

times. Moreover, in most models, the photons from the redshifted background are sufficient 
to fully reionise He II at redshift z ~ 3 without any additional contribution from quasars at 
lower redshifts and the temperature of the mean density intergalactic medium remains close 
to 10 4 K down to redshift z ~ 1. 
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1 INTRODUCTION 

The mean transmitted flux along the line of sight of the furthest 
quasa rs found by the Sloan Digital Sky Survey (e.g. jBecker et alJ 
1200 ll) shows that the ionisation fraction of the intergalactic medium 
(IGM) is rapidly decreasing approaching z ~ 6.3. Assuming, 
at that redshift, an approximately constant stellar emissivity of 
ionising radiation, the I GM would have been reionised at z le j 
~ 7 by ste llar sources 1 Gned irll2000iJ: iDioreovski eTai]l200ll: 
ISongaila & Cow ie 2002; Fan et al. 2003, and others). But a much 
earlier ionisation e poch is impli ed by the WMAP observations 
iBennett et alJl2003UKogut et alJl2003T) . If the IGM had a sudden 
transition from neutral to completely ionised, producing the op- 
tical depth to Thompson scattering r c = 0.17 ± 0.04 measured 



by WMAP, Zrei — 17 for the best fit ACDM cosmological model 
ISpereel et all2003h . 

In a companion paper (paper I. iRicotti & Ostrikerll2004 . wc 



have tried to explain these apparently contradictory results by con- 
sidering the contribution to reionisation due to zero-metallicity 
(Population III) stars. Since s uch stars are thought t o be mas- 



sive 1 


M, > 5 - 10 M n ) (e.e..lOmukai & Nishilll998tlAbel et al 


boor] 


Bromm et alj2002|:|Abel et alj2002ilNakamura & Umemura 


2002 


) and hot jTumlinson & Shullll200d). their efficiency for UV 



emission is larger than for Population II and Population I stars that 
have a Salpeter initial mass function (IMF). A time dependent ef- 
ficiency of UV emission, euv, produced by the transition from a 
top-heavy IMF at high redshift to a Salpeter IMF at low redshift, 
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could possibly explain the large r c measured by WMAP and would 
be consistent with the Sloan data at z £ 6.3 (Chiu et aL 20Q3). An 
upper limit of the UV emissivity, e™v x ~ 2 — 3 x 1CT 3 can be 
derived by assuming the maximum efficiency of energy produc- 
tion by thermonuclear reactions in massive stars. For a Salpeter 
IMF, euy ~ 2 — 4 x 10 ~ 4 , depending on the stellar metallicity 
(Tumlinson & Shull 2000). Therefore, assuming a constant (maxi- 
mal) escape fraction of ionising photons from the galaxies, (/esc) 
~ 1, a redshift dependent IMF would increase the UV emission 
efficiency from low to high redshift by about 10-20 times. 

Several authors have tried to reproduce the high optical depth 
measured by WMAP wit h the increa sed ionising photon emission 
by zero-metallicity stars. Cen 1 2003) has used a semianalytic cal- 
culation to indicate that the IGM ma y have experien ced two dis- 
tinct epochs of reionisation. lHaiman & HoldeJ {2003) have inves- 
tigated different reionisation histories showing that optical depths 
consistent with WM AP can be a chieved assuming very massive 
Population III stars. IWvi the & Loebl 12003) show that, if feedback 
regulates star formation in early low-mass galaxies as observed in 
present-day dwarfs, Population III s tars forming with a heavy IMF 
are required to match WMAP data. Ciardi et al. 1 2003), using nu- 
merical simulation, show that zero metallicity stars with a mod- 
erately heavy IMF (characteristic stellar mass of 5 M©) and (/esc) 
= 20% are able to generate an optical depth of 0.17, consistent with 
WMAP data. Sokasian et al. (2004), also using numerical simula- 
tions, showed that zero metallicity stars with a top-heavy IMF and 
(/esc) ~ 100% is needed to produce an optical depth consistent 
with WMAP data. 

In paper I we have found (consistent with the work of the other 
investigators mentioned above) that, assuming the most extreme 
properties for zero metallicity stars (top heavy IMF and (/ eS c) 
~ 100%), the maximum optical depth to Thomson scattering pro- 
duced by Population III stars is r c ~ 0.13, marginally consistent 
with WMAP data. In order to produce this large value of r c , mas- 
sive zero-metallicity (Population III) stars must be the dominant 
population until redshift z ~ 10. But we have noted that, if the ra- 
tio of metal atoms to ionising photons produced by Population III 
stars is normal, the metal enrichment of the ISM prevents metal- 
poor stars from being produced for long enough to ionise more 
than a small fraction of the IGM, z'.e.fhat the Population III phase 
is very rapidly self-limiting. In addition, if pair-instability SNe are 
important, the mechanical energy input by SN explosions produces 
strong outflows in galaxies with masses Mam ~ 10 9 M®, reducing 
their star formation and delaying reionisation. These two effects 
will lower the redshift of reionisation and r c to values inconsistent 
with WMAP but still consistent with the Sloan quasar data. Thus 
we have found it essentially impossible to reproduce the WMAP 
and Sloan results using stellar sources of UV ionisation. 

In summary, the reason for the partial disagreement with the 
conclusions of previous works on the ability of Population III stars 
to produce the large optical depth observed by WMAP, is motivated 
by our findings concerning the self-termination of Population III 
stars by the enhanced metal pollution and the negative feedback on 
star formation from the enhanced energy input by SN explosions, 
that inevitably follow from the assumption of a top-heavy IMF. In 
previous works these concerns have not been addressed. However 
in paper I we also note that, if most Population III stars collapse 
into black holes (BHs), the metal pollution problem is alleviated 
and the mechanical feedback from SN/hypernova explosions will 
be reduced as well. An interesting consequence of th is scenario is 
a cop ious production of BHs from the first stars l Schneider et al. 
2002). For this reason it is possible that the secondary radiation 



from accretion onto seed BHs might be a more important source of 
ionising radiation than the primary radiation from the same stars. 
In particular because, as we will show, accretion on seed BHs is 
not sensitive to the early termination of the Population III epoch by 
metal pollution. 

Motivated by these results, in this paper we study the partial 
ionisation of the IGM by an early X-ray background produced by 
accreting BHs, followed by more complete stellar reionisation by 
Population II stars at z rc i ~ 7. An hybrid model where UV from 
Population III stars and X-ray emission from accreting BHs are 
both important in producing a large optical depth is also viable if 
self-termination of Population III stars happens late (e.g., z ~ 10) 
and if we assume extreme parameters for (/esc) and for the global 
Population III stars UV emissivity. But, given the difficulties in 
modelling properly the transition from Population III to Popula- 
tion II stars, the predicted importance of UV from Population III 
stars is very uncertain (almost a free parameter at this point). As an 
example, in the present paper we show a model where the transition 
from Population III to Population II stars is at redshift z ~ 15. In 
this case we show that UV radiation from Population III stars has 
negligible effect on the Thompson optical depth. 

For a fixed emissivity of ionising radiation, X-rays are less ef- 
ficient in reionising the IGM, because most of their energy goes 
into heat instead of ionisation. The inclusion of secondary elec- 
trons can boost the ionisation efficiency by about a factor of ten if 
the electron fraction of the IGM is (x c ) < 10 %, but roughly 1/3 of 
the primary electron energy is always converted into heating of the 
IGM JShull & van Steenbergll 19851) . It will be shown that the par- 
tial ionisation of the IGM must begin at redshifts z ^ 15 in order 
to have an important effect on r c . This can only happen if the X-ray 
sources (i.e. , accreting BHs) form in small-mass haloes with masses 
Mj m ^ 10 s Mo . Previou s work on the formation of the first galax- 
ies l R icotti et all20023bl) has used cosmological simulations with 
radiative transfer to show that star formation in the first small-mass 
galaxies is reduced by feedback effects but is not suppressed. It is 
therefore plausible that a substantial production of seed BHs from 
the first stars took place in the first small-mass galaxies. 

Published studies have investigated the effects of X- rays 
on ea r ly galaxy format i on usi ng a semianalytic approach lOhl 
l200lt I Venkatesan^ydJ 1200 ll) and cosmological simulations 
I Mac haceket all 2003). But since those works were completed be- 
fore WMAP results, they did not focus on scenarios that could 
produce r c ~ 0.17. The study presented in this paper (pa- 
per Ha) is carried out using a semianalyti c code bas e d on p rin- 
ciples similar to those adopted by IChiu & OstrikeJ l2000l) . In 
iRicotti. Ostriker. & Gn edin 1 200% (paper lib) we present the re- 
sults of cosmological hydrodynamic simulations that include ra- 
diative transfer. We discuss the observational signatures of X-ray 
preionisation compared to stellar reionisation, including calcula- 
tions of the expected amplitude of secondary anisotropies of the 
CMB and the redshifted 21cm signal in absorption and emission 
against the CMB. The semianalytic models presented in this paper 
allow us to explore a larger parameter space and help us to inter- 
pret the numerical result of the full cosmological simulations. This 
new code, presented in Appendix^] calculates the evolution of the 
filling factor of H II regions and the ionisation and thermal history 
of the IGM outside the H II regions, produced by an X-ray back- 
ground. The semianalytic code has been tested and calibrated using 
the results of the cosmological simulations (cf, paper I). Note that 
some cosmological simulations presented here include the effects 
of SN explosions, using the recipe dis cussed inlGnedinl 1199 8). in 
contrast to earlier papers in this series l Ricotti etal. 2002jbh. 
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This paper is organised as follows. In §[2] we present qualita- 
tive arguments showing why and for which models X-ray preion- 
isation can be more effective in increasing r e than reionisation by 
stellar sources. In §|5|we estimate the accretion rate onto seed BHs 
and we compute the integrated X-ray energy required to explain 
the large r c measured by WMAP. In §[4] we summarise the char- 
acteristics of the semianalytic code for reionisation and we show 
its results. In § [5] we conclude this work with a summary and a 
discussion of some of the observable effects of early black hole X- 
ray heating. A full description of the semianalytic code is given in 
Appendix IaI 



2 X-RAY PREIONISATION: RATIONALE 

The essential physics is easily understood. Ultraviolet (UV) pho- 
tons from stars at wavelengths only moderately shorter than 
912 A (hv = 13.6 eV) have a very short mean free path in neutral 
hydrogen. As a consequence they produce Stromgren spheres sur- 
rounding the UV sources in high density regions. In such regions 
recombination rates are high and many photons are required to keep 
one hydrogen atom ionised for a Hubble time. The much larger 
volumes between the Stromgren spheres remain sensibly neutral. 
However, much harder photons {e.g., 1 keV) can uniformly fill 
space and thus lead to fractional ionisation which is higher in the 
low-density regions than in the high-density regions containing the 
sources. But, as noted earlier, the fractional ionisation is low but 
the volume is large since the volume fraction in virialised haloes at 
z = 15 is only a few per cent. To illustrate this let us show how a 
simple ad hoc ionisation history affects r c . In Fig.Q we show the 
IGM optical depth to Thompson scattering, r c , in a model where a 
complete reionisation at redshift z rc i is preceded by a partial ioni- 
sation, with constant ionisation fraction (x e ), starting at z = zpbh- 
The long-dashed, dashed, dotted and solid lines in Fig.Qshow r c as 
a function of zpbh assuming complete reionisation at z rc i = 7 and 
preionisation electron fraction (x c ) = 0.7, 0.5, 0.3 and 0.2, respec- 
tively. We see that two conditions need to be met in order to have 
t c ~ 0.17: (i) preionisation must start early, at zpbh ~ 20 — 40, 
and (ii) the IGM preionisation fraction must be (x c ) £ 20%. 

Condition (i) is met only if the X-ray sources can form in 
haloes with masses A/dm ~ 10 8 M© . The two peculiarities of these 
small-mass haloes are that they rely on molecular hydrogen cool- 
ing to form stars, and their dark matter gravitational potential is 
too shallow to retain gas hotter than T w 2 x 10 4 K. For these 
reasons small-mass galaxies are subject to ra diative f eedback ef- 
fects that might suppr ess their formation (e.g.. | Haiman et alll997t 
Mach acek et alj F2003. but see lRicotti et alj|200ll who find a new 
positive feedback effe ct). Us i ng cosmo logical simulations with ra- 
diative transfer Ricotti et al. 1 2002a b) found that the main nega- 
tive feedback thought to suppress the formation of the first galaxies 
(the H2 photodissociating background) is not the dominant effect. 
Instead, feedback by ionising radiation plays the dominant role, 
inducing a bursting star formation mode in small-mass galaxies. 
The feedback mechanism is complex because it is composed of 
several processes, each with different relevance depending on the 
parameters of the simulation (mainly the ionising escape fraction, 
(/esc), and the IMF). In brief, photoevaporation and H2 forma- 
tion/photodissociation produce a bursting mode of star formation 
that is catalyst for molecular hydrogen re-formation inside relic H II 
regions and in the precursors of cosmological Stromgren spheres. 
Star formation in small-mass galaxies is self-regulated on a cosmo- 
logical scale and it is reduced by radiative feedback but is not sup- 
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Figure 1. Optical depth to Thomson scattering, r e , for a toy model where 
a complete reionisation at redshift z rc i = 7 is preceded by constant, par- 
tial ionisation (x c ), starting at redshift z = 2pbh- The lines from bot- 
tom to the top refer to mean (mass weighted) ionisation fractions (x e ) 
M = 0.2,0.3,0.5, and 0.7. 



pressed. Even if the volume fraction filled with ionised Stromgren 
spheres never exceeds a few per cent, these small-mass galaxies 
form enough Population III stars and therefore seed black holes 
(see §|3J to produce, via accretion, a substantial X-ray background 
at high redshift. 

Condition (ii) can be met if the X-ray energy input is suffi- 
ciently large. Note that ionisation by secondary electrons becomes 
increasingly inefficient when (x e ) Si 10—20% since the secondary 
electrons loose an increasing fraction of their energy to heat and a 
decreasing fraction to further ionisation as the neutral fraction de- 
clines. We can estimate the fraction of baryons, uj ac , that we need to 
accrete by counting the number of ionisations per hydrogen atom, 
iVxph — (x e ), needed to partially ionise the IGM to a level (x e ) 
~ 50% (the fractional ionisation of the IGM must be (x a ) £ 20% 
in order to produce the large optical depth to Thomson scattering 
measured by WMAP, cf. Fig.QJ. If we neglect recombinations we 
have 

iVxph = 7.3 x 10 7 ( ) cj rc ~ (x c ) -0.5, (1) 

\ AXscc / 

where ex is the energy emitted in the X-ray band per hydrogen 
atom accreted on the BHs and Xxscc is a pure number defined 
as Axscc = X%/f sec ~ 20, where / scc ~ 5 — 10 takes into 
account the extra ionisation by secondary electrons (we have esti- 
mated this number switching off secondary ionisations in our semi- 
analytic code, presented in the appendix), and Xx = hv/(13.6eV) 
is the mean energy of X-ray photons in units of the H I ionisation 
energy for the assumed X-ray spectrum. For accretion around a BH 
in which the gas is highly ionised and electron scattering is the 
most important opacity, the isotropic luminosity produced by the 
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accretion cannot exceed the Eddington luminosity 

_ 4nGM B nmp 38 _ x / M B h\ 

L E d = = (1-3 x 10 erg s ) — — , (2) 

a T V M Q J 

where ot is the Thomson scattering cross section. Given the ac- 
cretion efficiency e = Z//M ac c 2 ~ 0.15, typical for quasars, 
the mass accreted at the Eddington rate is MEd = Z/Bd/ec 2 
and the timescale of the accretion is the Eddington time t^d = 
M BH /M E d ~ 10 s yrs. 

Observations of active galactic nuclei (AGN) at low-redshift 
show that, when the central BH is accreting, it does so at near the 
Eddington limit and, when it is quiet, the energy output is negli- 
gible. In §|3 we will introduce a parameter, /duty, to model the 
observed duty cycle of AGN. We define /duty as the fraction of 
time when the BH is accreting at the Eddington limit. Assuming a 
spectral energy distribution (SED) typical for AGNs (cf., Fig. [4}, 
a fraction = Lx/^Ed ~ 10% of the energy output from the 
accreting BH is emitted in the X-ray band, with a resulting effi- 
ciency of X-ray production ex = (3e. Finally, using equation Q, 
we can estimate the fraction of baryons that needs to be accreted to 
partially reionise the IGM, 

'0.02\ f X Xse , 
ex J V 20 

As a comparison we also estimate the baryon fraction needed to be 
converted into stars to reionise the IGM. The maximum radiative 
efficiency from thermonuclear processes is £uv x — 2 — 3 x 10 -3 . 
Since UV photons have a short mean free path they first ionise, 
and keep ionised, the high density regions around each source be- 
fore escaping into the IGM. In this case, then we cannot neglect 
recombinations. The number of UV photons per baryon needed 
to reionise the IGM is Vuvph — Chu 10, where Chu = 
(n 2 ) / (^c) 2 is the effective clumping factor of the ionised gas. Note 
that the effective clumping, Ch □ , is the time averaged clumping of 
the IGM from the time the first source turns on to the time of reion- 
isation. When the volume filling factor of the Stromgren spheres 
is small the clumping factor of the ionised gas is very large, but 
equals the clumping of the gas, C, at the redshift of reionisation 
(see Fig. 4 in paper I). In addition, the escape fraction of ionised 
photons, (/esc), and the value of Chu are related to each other and 
depend on the resolution of the hydrodynamic simulation or on the 
definition of (/ CS c) in semianalytic models. Keeping in mind these 
caveats, the fraction of baryons that needs to be converted into Pop- 
ulation III stars to reionise the IGM is approximately 



6.8 x 10" 



(3) 



lu* ~ 4x10" 



Chii 

To" 



0.5 



2 x 10" 
euv 



A't 



(4) 



(/esc) J V £ UV J V 3 

where Auv = hu/(13.6eV) is the mean energy of UV photons 
in units of the H I ionisation energy for a Population III spectrum. 
By comparing equation|4|to equation|3it is evident that the baryon 
fraction that needs to be accreted onto BHs to partially reionise 
the IGM to (x c ) — 0.5 is about 60 times smaller than the fraction 
of baryons that needs to be converted into massive Population III 
stars. This shows that X-rays from accreting BHs are potentially 
more efficient than UV from stars in producing a large r c . 

Observations in the soft X-ray bands of Lyman-break galaxies 
(which do not contain AGNs) at redshift z ~ 3 iNandra et a l. 2002) 
have shown that their soft X-ray emission is proportional to the UV 
emission: Lx/iuv ~ 0.01. This emission, attributed to X-ray bi- 
naries and SN remnants, is too low to influence preionisation. In 
addition, about 3% of Lyman-break galaxies contain central AGNs, 
accreting at roughly the Eddington limit l Steidel et al. 2002J). X- 
ray emission powered by accretion on primordial seed BHs could 



thus be important. There are two scenarios for the formation of 
supermassive black holes (SMBHs) in the bulges of today's galax- 
ies: monolithic collapse, or mergers and accretion on seed BHs. 
Only the mass growth due to accretion produces X-rays available 
to ionise the IGM. Assuming that most of the mass in SMBHs ob- 
served in the bulges of galaxies at z — has been accreted at high 
redshift, the upper limit for the accreted mass is 



WBH.max 



/ MsMBH \ 
V -^bulgc / 



^buigc = (1 ± 0.7) x 10" 



(5) 



where we used MsMBit/Mbu 



i ge = (1-5 ± 0.3) x 10" 3 (e.g., 
lKormendv"& Richston e ^^9^|G£bhard^^lj2^ oT) and uJh„i^ = 
(6.5 ± 3)% frersic & Saluccll992l:lFukugita et alll998t) . But an- 
other scenario is possible in which BH seeds grow through accre- 
tion at high redshift to a level higher than currently observed, but 
a significant fraction of them is expelled from the galaxies or does 
not end up in the observed population of SMBHs. In this case a 
fraction of the halo DM would be made of primordial BHs. The re- 
cent discovery of several ultra-luminous X-ray sources (ULX) has 
been interpreted as evidence for intermediate-mass BHs (about 100 
Mq). These objects could also be relics of primordial BHs. 

In summary, X-rays may be more effective in increasing the 
optical depth to electron Thompson scattering r c with respect to 
UV photons from stellar sources for the following reasons: 

• X-rays escape unabsorbed from the ISM of galaxies while for 
UV radiation (/esc) < 1. 

• X-rays ionise preferentially the low-density regions where the 
clumping factor is Chu < 1. Consequently, the recombination rate 
is reduced. 

• X-rays only partially ionise the IGM: (x c ) ~ 20 - 70%. This 
reduces the recombination rate by a factor (x e ). 

• The emission efficiency of radiation due to gravitational ac- 
cretion can be as high as e ~ 0.15 (accretion at the Eddington 
limit). This efficiency is about 100 times the maximum efficiency 
produced by thermonuclear reactions in stars ejjv 1 — 2 x 10~ 3 . 

• The accretion on seed BHs can continue even if the epoch of 
Population III star domination is very short due to the metal pollu- 
tion of the star forming gas. 

On the other hand, X-rays may be inefficient with respect to UV 
radiation in ionising the IGM for the following reasons: 

• X-rays have to be emitted at early times to be effective. This 
can only happen if a substantial density of seed BHs is formed in 
small-mass galaxies, at z £ 20. 

• X-rays are less efficient in ionising than UV photons. But, if 
the electron fraction in the IGM is (x e ) ~ 10 — 20 %, secondary 
electrons are an important additional source of ionisation. Never- 
theless, as we will show in § 14.21 hard-UV photons and the red- 
shifted X-ray background can enhance the IGM electron fraction 
to levels well above (x e ) ~ 10 — 20%. 

• X-rays are efficient in heating the IGM and therefore they in- 
crease the IGM Jeans mass. This will have a negative feedback on 
the formation of the smaller mass galaxies and consequently on the 
formation of seed BHs and their ability to accrete gas. 



3 MASS ACCRETION HISTORY ON SEED BHS 

The formation rate of seed black holes from star formation is very 
uncertain. The calculations of the mass of the remnant and metal 
yields are strongly dependent on the energy of the SN explosion, 
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Figure 2. (left) The top panel shows the black hole accretion rate, p ac = po<i> ac (where po is the cosmic baryon density at z = and o; ac is given in 
equation Q), as a function of redshift. The bottom panel shows the baryon fraction in BHs, ubh, as a function of redshift. (right) The top panel shows the 
star formation rate, p* = pout*, as a function of redshift. The bottom panel shows the baryon fraction in stars as a function of redshift. The solid, dotted and 
dashed lines refer models M3, M4, and M5 in table[2] respectively. The dot-dashed line in the top panel of the left figure is an upper limit for the accretion rate 
given by equation 1111 . 



£51 = E/10 erg, on the mass shell at which the explosion takes 
place and the star metallicity. These numbers are uncertain and cal- 
culations ( Wooslev & Weavea l 19951) have been done for a grid of 
models. In tableQwe list the fraction, /bh, of stellar mass expected 
to end in seed BHs for the different explosion models, gas metallic- 
ides and stellar IMF. The data is from lWooslev & WeaveJ ll995l) : 
Model A has E51 ~ 1; Model B has E51 ~ 2; and Model C has 
£51 ~ 3. The values listed in the table are the weighted means 
of /bh for stellar masses between 12 Mq < M» < 40 Mq with 
weight M a . Note that for a Salpeter IMF (a = -1.35) with stellar 
masses between 1 Mq < M* < 100 Mq, the mass in stars with 
M» > 8 M is 20 % of the total. 

Very massive stars (VMS) in the mass range 140 Mq < 
M* < 260 M are believed to explode completely without leav- 
ing any remnant due to an instability produced by annihilation 
of electrons and positrons (pair-instability SNe) and hence have 
/ BH = 0. Stars more massive than M* = 260 Mq are believed 
to collapse directly into BHs, without exploding as SNe and hence 
have / BH = 100% <Heeer & Wooslevll2002l) . The existence of a 
mixture of masses for supermassive Population III stars allows for 
all possible values of /bh from zero to unity. 

The mass growth of seed BHs can be calculated if we know the 
accretion rate w ac of BHs, the formation rate of seed BHs, w see d, 
and the ejection rate of BHs from galaxy haloes, w e j • To calculate 
the mass growth of BHs we solve the following equation: 

WBH = <^ac — <^cj + <jj S ocd , (6) 

where ojbh, k>ac, k> c j and ui scc d are expressed as a fraction of the 
baryon content of the universe. The accretion rate is proportional 
to wbh, independently of the mass function of BHs, 

<jJbh 



Table 1. Fraction of stellar mass that forms seed BHs 
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7 


1 


10" 4 


10 


g 


8 


-1.35 


10" 4 


13 


8 


7 





10~ 4 


16 


8 


7 


1 


1 


11 


10 


8 


-1.35 


1 


13 


to 


8 





1 


15 


to 


8 


1 



where the accretion time, i ac = tEd//duty. is expressed as the 
fraction of time, /duty, that the BHs are accreting at the Eddington 
rate. If we assume that the BH accretes at the Eddington limit for 
a time interval t on and stops accreting for a time interval t a «, the 
duty cycle parameter is /duty = ton /(toff + ion) < 1. If we neglect 
BH ejection we have 



, ^BH , 

WBH = /duty— 1 /BHW,, 

lEd 



(8) 



where we have assumed that the formation rate of seed BHs is a 
fraction, /bh, of the star formation rate u ECC d = /bh<^« (see ta- 
ble|2j- Integrating equation JS} we get 



where, 



/duty 



^Ed 



(7) 



lobh (t) = f di'/BH^, (*') exp [F(t) -F(t% 
Jo 



F(t) = — / fdutydt'. 
tEd Jo 



(9) 
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We parameter! se the duty cycle as 



Table 2. Parameter of the semianalytic models 



/duty (z) 



l + Z 
ZBH 



(10) 



and imposing the condition 0.001 ^ /duty ^ 1 to be consis- 
tent with the obser ved fraction of AGN /duty ^ 3% at z ~ 3 
ISteidel et all2002l) . and f duty ~ 0.001 at z = 0. 

In this work we will not attempt to investigate in detail the 
physical processes and feedback effects that could produce the duty 
cycle given in equation <10l . Here we simply point out the exis- 
tence of, at least, one simple physical model that can reproduce 
the family of equations given in J10I . In this model the time i ff 
during which the BH is quiescent is the Compton cooling time 
tcomp = (224 Myr)[(z + 1)/10]~ 4 . This is the time that it takes 
for the gas in the proximity of the BH to cool and be accreted after 
the temporary halt of the accretion produced by gas heating from 
the previous BH activity. If we further assume that the duration of 
the burst is a fraction, rj, of the Eddington time (t on = ^Ed) we 
can reproduce /duty = 1/(1 + toff Aon) with the form given in 
equation 1101 . where the value of zbh is determined by the value 
of rj. For reasonable values of rj in the range 1 — 100% it is possible 
to describe models spanning from the early (i.e., adopting rj ~ 1%) 
to the late preionisation (i.e., adopting rj ~ 100%). On the basis of 
observations of AGNs at redshift z £ 3, the best estimate of t on is 
~ 10 Myr (e.g.,|Haehnelt et al. 1998). This suggests that the inter- 
mediate preionisation models (that adopt rj ~ 10%) are prefered. 

For the purpose of this paper we do not need to focus on the 
functional form of the duty cycle at low redshifts (z £ 3) provided 
that we can avoid to violate observational constraints on the duty 
cycle and BH accretion rate. In order to get a realistic accretion rate 
at lower redshift in some models, it might be necessary to adopt a 
non zero ejection rate of BHs from the host galaxies [i.e. , w e j ^ in 
equation J6j]. Also this assumption has physical motivations (e.g., 
Mad au et alj2004 . 

In §|4|we use a semianalytic code to follow the chemical and 
thermal history of the IGM in models where accretion on seed BHs 
partially ionise the universe. At redshift z ~ 6.3 stellar sources 
reionise the IGM to a level consistent with Sloan quasar observa- 
tions. We explore three cases that differ for the accretion history 
onto seed BHs. The accretion history, parametrised using the equa- 
tions derived in this section, can happen early (z ~ 25), late (z ~ 
15) or at intermediate redshifts (z ~ 20). If we wished to allow for 
a lag, At, between the formation of an individual BH and the com- 
mencement of accretion on it (cf., Mada u et all2004IWhalen et alJ 
120031) we would simply replace w*(t') with u)*(t' — At) in equa- 
tion j9j- In table |2| we list the parameters adopted for the semian- 
alytic simulations whose results will be shown in §0 The values 
for /bh allow for a range of mass functions from Salpeter to top- 
heavy. We also set /bh = in all models after z = 15 since the 
formation of Population III stars that produce the seed BHs is prob- 
ably self-terminated by metal enrichment well before reionisation 
at z rc i = 7. But this assumption does not have any effect on the 
BH accretion history because the total mass in BHs becomes domi- 
nated by the accretion on e-folding time scales t ~ t ac , making the 
contribution by further seed BH formation negligible after the ini- 
tial phase (before the global accretion rate reaches the maximum). 

The BH accretion rate as a function of redshift, p ac (z) = 
Po^ac, where po is the cosmic baryon density at z — and o; ac is 
given in equation {7), is shown in the top panel of Fig.|5Jleft). The 
baryon fraction in BHs, ljbh, as a function of redshift is shown 
in the bottom panel of Fig. [2[left). The solid, dotted and dashed 



Model 


/bh (%) 


4> 


ZPBH (*on Myr) 


n s 


M3 


0.2 


12 


15 (~ 100) 


1.00 


M3b 


0.1 


8 


15 (~ 100) 


1.04 


M4 


2 


8 


20 (~ 10) 


1.00 


M4b 


2 


8 


20 (~ 10) 


1.00 


M5 


20 


12 


25 (~ 1) 


1.00 



■ZPopHI 



15 



Meaning of the parameters: /bh is the fraction of stellar mass that col- 
lapses into BHs; <j> and zpbh are the parameters in equation llOl that de- 
termines the time evolution of the duty cycle /duty; n s is the spectral 
index of the spectrum of primordial density perturbations; zp op in is the 
redshift at which UV radiation from Population III stars (that has emissiv- 
ity 1 1 times larger than the UV from Population II stars) begins to fade: 
the function for the UV emissivity is given by £uv( 2 )/ £ uv( 2 = 0) = 
10{arctan [2(z — zp O plll)]/ 7r +0.5} + l, where the emissivity from Pop- 
ulation II stars (that is used in all the other models) is euv( 2 = 0) = 
4.8 X 10 -5 . 



lines refer models M3, M4, and M5 in table|2| respectively. As al- 
ready noted, all the models have about the same ojbh ~ 10" 4 at 
z = (i.e., the same total energy output from the sources) and 
differ for the characteristic epoch of accretion. The BH accretion 
rate is always negligible at the redshift, z rc i, of IGM reionisation 
by stellar sources. The dot-dashed line shown in the top panel of 
the Fig.|2jleft) is an upper limit for the accretion rate. This limit on 
the accretion rate is given by equation Jl II . that will be discussed 
in5l4"2l 

In Fig.|3right) we show the star formation rate (top panel) and 
the baryon fraction in stars, w» (bottom panel), as a function of red- 
shift for the same models (cf., table[2j. The p oints with errorbars i n 
the top panel show the observed SFR from lLanzetta et aljEo53) . 
The baryon fraction in stars at redshift z = is u) t = 20% in a ll the 
models, consistent with observations 1 Pers ic & Saluccil 19921) . The 
feedback due to the X-ray heating of the IGM suppresses galaxy 
formation in haloes with masses smaller than the IGM filtering 
mass, producing the substantial reduction of the global SFR ob- 
served at redshifts around z = 15 in the top panel of Fig. fright). 

Finally, in Fig. [5] we show the mass fraction accreted onto 
seed BHs with respect to the star fraction in the universe, p a c/p* = 
aj ac /aj* (top panel), and the adopted duty cycle (bottom panel). 
Again, the dashed, dotted and solid lines refer to models M3, M4 
and M5 in table |2| As shown by the labels next to each line, the 
early accretion model has /bh = 20%, therefore it requires that 
most of the mass of Population III stars implodes into BHs. But 
for the other models only a fraction /bh 2% of the stars needs 
to implode into BHs. This condition does not require a top-heavy 
IMF. It is interesting to note (but not surprising) that in order to 
have a substantial mass accretion onto seed BHs at early times a 
large fraction of seed BHs, /bh, from Population III stars is re- 
quired. This is because the accreted mass increases exponentially 
with e-folding time scale igd — 10 8 yrs, which is about the Hub- 
ble time at z = 30. To conclude we note that if seed BHs are not 
able to accrete initially at the Eddington rate (see bottom panel of 
Fig.[3]that shows the adopted duty cycle), a larger fraction of seed 
BHs, /bh, is necessary in order to get the same global accretion 
rate as in our models. This is clearly possible to do for the interme- 
diate and late accretion models, but would need extremely high and 
unrealistic values of /bh in the early accretion model M3. 
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Figure 3. The top panel shows the mass fraction accreted onto seed BHs 
with respect to the star fraction in the universe. The bottom panel shows 
the duty cycle, /d u ty> assumed in each model. The solid, dotted and dashed 
lines refer models M3, M4, and M5 in table|2] respectively. 



3.1 Spectral Energy Distribution 

In Fig. [4] we show the spectral energy distribution (SED) of stars 
and quasars adopted in our simulations. For the stars we adopt 
a SE D calculated for zero-metallicity stars fTumlinson & ShuTH 
2000) and for quasars w e use a spectrum simi lar (but not identical) 
to the one calculated bv lSazonov et alj 120041) . Their spectrum is a 
template, inferred on the basis of observations, representing the av- 
erage quasar in the universe and determined by considerations both 
of the (X, 7) ray background and the spectra of individual objects. 
The quasar spectrum has an X-ray bump produced by absorption 
of UV photons by intervening gas and an IR excess produced by 
the reprocessed UV photons. Using this SED we can distinguish 
between reionisation by UV photons produced mainly by the stars 
and X-rays produced by quasars. The energy band that contributes 
the most to reionisation, also amplified by the effect of secondary 
electrons, is 0.5 keV < h v v < 5 keV. But as we will point out in 
paper lib, the X-ray background photons that are redshifted to the 
extreme ultraviolet (EUV) bands can be the dominant source for 
preionisation in the voids. 

The typical expected mass of the accreting BHs in our model 
is M B h £ 100 - 1000 Mq at the high redshifts of interest, while 
the adopted SED is inferred from the accretion onto BHs with 
masses of M BH £ 10 5 - 10 6 M Q . Generally speaking the SED 
of accreting compact objects has two components that contribute 
comparably to the total energy output: (i) a power law compo- 
nent from the hot corona around the BH and (ii) a thermal com- 
ponent from the inner part of the hot accretion disk. The power 
law component is independent of the mass of the accreting ob- 
ject (from X-ray binaries to AGN) while the disk component has 
a temperature that depends mildly on the mass of the BH as 
fcT disk « 1.2 keV(M BH /10 M Q )- 1/4 . The observed spectra of 
ULX sources tentatively show a subdominant soft component fit- 
ted by multicolour disks blackbody with fcTdisk ~ 0.15 keV, 




h p i/[Kev] 

Figure 4. Spectral energy distribution of Population III stars (dashed line) 
and quasars (solid line). UV radiation is mainly emitted by stellar sources 
and X-rays by quasars. We assume that most UV radiation from quasars 
does not escape to the IGM but it is absorbed by an obscuring torus or in 
the ISM. In the adopted quasar SED the thermal component, that depends 
mildly on the BH mass, is negligible with respect to the power low compo- 
nent that is mass independent. 



which has been interp reted as a tent ative evidence of intermediate 
mass black holes 1 Miller et al. 2003). The non-thermal component, 
which dominates the total energy output, is a power law similar to 
the one adopted in the present pap er. Note that the spectrum that we 
have assumed is also similar to the lSazonov et alJl2004l) "absorbed 
quasar spectrum" (Fig.[4j were the UV bump produced by the ther- 
mal component in supermassive BHs is absorbed and in fact has 
no UV at all. The spectrum that we adopt is therefore dominated 
by the power law component (absorbed in the UV) that is roughtly 
independent of the assumed mass of the accreting BH. 



4 SEMIANALYTIC MODELS 

In this section we show the results found using a semianalytic 
model for reionisation. The advantage of the semianalytic approach 
with respect to the cosmological simulations presented in paper lib 
is that we can explore a larger parameter space and the results 
are not affected by the resolution of the simulations. We use this 
approach to study the dependence of the Thomson optical depth 
on cosmological parameters and accretion histories of seed BHs. 
The results presented in this section also helped to derive a phys- 
ically motivated, time dependent X-ray emissivity, to be used in 
the more computationally expensive numerical simulations. In ad- 
dition, the semianalytic models are an aid for the interpretation 
of the results of cosmological simulations, made complicated by 
the interplay of a larger number of physical processes. If not oth- 
erwise specified, we adopt the concordance ACDM cosmologi- 
cal model with h = 0.72, fi m = 0.27, Q h = 0.046, n s = 1 
and (78 = 0.91, that is consistent with the WMAP measurements 
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Figure 5. (left) The top panel shows the ionised hydrogen fraction and the volume filling factor of H II regions (dot-dashed line) as a function of redshift for 
model M3. The bottom panel shows the IGM temperature as a function of redshift for model M3. The dotted, solid and dashed lines refer to gas over-densities 
6 = 0.1, 1, 10 in the IGM outside Hll regions, respectively. After overlap, when the ionising background is uniform, the lines refer to the ionised IGM. 
(centre) Same as in the left figure but for model M4 (thick lines) and M4b (thin lines), (right) Same as in the left figure but for model M5. 



(h = 0.71 t^t ttmh 2 = 0.135±^, Q ^h 2 = 0.0224±0.0 009, 
n s = 0.93 ± 0.03 and cr 8 = 0.84 ± 0.04. ISpergel et alj2003t) . 



4.1 The code 

We implemented a semianalytic model to study IGM reionisation, 
chemical evolution and re-heating. The filling factor of H II regions 
is calculated following the method in IChiu & Ostrikeil feOOd) . But 
we also include X-ray ionisation of the IGM before overlap of the 
Hll regions. A detailed description of the code is given in Ap- 
pendix^] Here we summarise the main processes included in the 
code at this stage: 

• We calculate the mass function of DM haloes and the forma- 
tion rate of haloes using the extended Press-Schechter formalism. 
The star formation rate (SFR) is proportional to the formation rate 
of haloes. In the integral for the global SFR we take into account 
cooling and dynamical biases that depend on the mass of the haloes. 

• We include thermal feedback on star formation: the minimum 
mass of star-forming galaxies at a given redshift is determined by 
the IGM filtering mass. 

• The emissivity of accreting BHs is calculated as in §|5| 

• We assume a two-phase IGM. One phase is the fully ionised 
gas inside the H II regions and the other is the partially-ionised or 
neutral gas outside the H II regions. 

• Radiative transfer is solved by splitting the spectrum into 
H I optically thick radiation and background radiation. Given the 
clumping factor inside H II regions we follow the evolution of the 
filling factor, temperature and chemistry of the ionised gas. Out- 
side the H II regions the temperature and element abundances are 
calculated as a function of the overdensity. 

• We calculate the time-dependent chemical network for 8 ions 
(H, H + , H~, H 2 , H+, He, He+, He ++ ) and thermal evolution in- 
side and outside Hll regions. The following heating and cooling 
processes are included: collisional and photo-ionisation, ionisation 
by secondary electrons, H, He and H2 cooling, Compton cooling- 
heating and the cosmological expansion term. We neglect cooling 
from metal contaminatio n in small-mass galaxies. The rates are 
from lRicotti et alJ feOOll) . 



4.2 Results 

We have run several models to study the sensitivity of the results to 
the accretion history on seed BHs and to cosmological parameters. 
In table [2] we list the parameters of the models shown in the next 
figures. In all the models the integrated energy emitted by accreting 
BHs is the same and stellar reionisation happens at 2 ~ 7. We nor- 
malise the BH accretion rate assuming that the total mass in BHs at 
2 — is WBH,max = 10~ 4 , the mass of SMBHs in the galactic cen- 
tres. This assumption is rather conservative since it is likely that a 
fraction of seed BHs is not incorporated into the SMBH population 
at z = 0. During galaxy mergers there is a non-negligible probabil- 
ity that one of the two merging BHs sitting in the bulge is expelled 
from the galaxy (cf., Mada uet al 12004)) . This probability is larger 
for unequal mass BHs, and happens in the final phase of the merger 
during the phase where gravitational radiation is emitted. 

In summary the constraints on the semianalytic models are the 
following: 

• The integrated energy output from accreting BHs is limited 
by assuming that the total mass in BHs at 2 = equals the total 
mass of SMBHs in the observed in the centres of galaxies. This 
requireme nt is a very conservative one, since theoretical calcula- 
tions (e.g. jMadau & ReeJ20 01) show that most intermediate mass 
BHs formed at high redshift will not merge into the SMBHs in the 
galactic centres. 

• Seed BHs accrete initially at the Eddington rate. The duty cy- 
cle in the late preionisation model (M3) is consistent with observa- 
tions of QSO at z ~ 3 that show a duty cycle of 3% and with QSO 
at z ~ that show a duty cycle of 0.1% 

• The star formation efficiency is constrained by fitting t he ob- 
served star formation rate (Madau plot) at z < 6 ( Lanze tta et alJ 
2002). The integrated baryon fraction in stars at z — is consistent 
with the one estimated by Fuku gita et alJ ll998l) . 

• The effective UV emissivity from stellar sources is con- 
strained to reproduce the transmitted flux observed in the Sloan 
quasars at z < 6.2. Assuming a Salpeter IMF the adopted escape 
fraction is (/ CS c} = 30 %. 

• The specific intensity of the ionising background at 2 < 6 is 
consistent with observations. 
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The typical ionisation and thermal history of the X-ray preion- 
isation scenario is illustrated in Fig.[3]for hydrogen and Fig.[6]for 
helium. In Fig.[5]we show the hydrogen ionisation (top panel) and 
the IGM thermal history (bottom panel) for model M3 (figure on 
the left), model M4 and M4b (figure in the centre) and for model 
M5 (figure on the right). In model M3, BH accretion takes place 
mostly at z > 15, at z > 20 in model M4, and at z > 25 in model 
M5. The dotted, solid and dashed lines refer to gas over-densities 
5 — 0.1, 1, and 10, respectively. Before reionisation at z rc i ~ 7, 
the lines refer to the gas outside the Hll regions. After overlap, 
when the ionising background is uniform, the lines refer to the fully 
ionised IGM. The dot-dashed lines in the top panels show the vol- 
ume filling factor of the H II regions as a function of redshift. In the 
figure at the centre the thin lines show a model (M4b) where the UV 
emissivity is about 10 times larger at z > 15 than at z < 15, to ac- 
count for the contribution of a top-heavy Population III star popula- 
tion. The volume filling factor of the H II regions in this model dif- 
fers from the same model where the Population III contribution is 
neglected (M4). The optical depth to Thompson scattering in model 
M4 and M4b has about the same value. 

A s even the earliest investigations showed ( Ostri ker & Gnedirl 
Il996l) . reheating significantly precedes reionisation (weather or not 
X-ray ionisation is included). In any X-ray preionisation model, by 
the time the ionisation fraction in the IGM reaches (x c ) ~ 10%, the 
IGM has a temperature T w 10 4 K (cf., Fig.[3J. This happens in all 
the models independently of the starting redshift of preionisation 
because the ionising efficiency of secondary electrons is reduced 
when (x c ) ^ 10% and the energy of X-ray photons is mostly con- 
verted into heat. 

By inspecting the dotted, solid, and dashed lines in the bot- 
tom panels of Fig. [5] that show the gas temperature at overdensi- 
ties S — 0.1, 1, and 10, respectively, we identify three epochs in 
the reheating history of the IGM. Initially, before complete reion- 
isation by UV radiation, the temperature of the IGM outside the 
cosmological Hll regions is almost isothermal (Tigm ~ 10 4 K) 
with overdense regions slightly cooler than underdense regions 1 . 
At reionisation the IGM becomes isothermal. After reionisation the 
IGM temperature dependence on the density follows a tight re- 
lationship, often referred to as "the effective equation of state of 
the I GM", Tigm = Tigm .o (1 + 5) (7_1) , where S is the overden- 
sity (Hui & Gnedin 1997). The parameter (7 — 1) is zero (i.e., the 
IGM is isothermal) at reionisation and (7 — 1) > 1 (i.e., over- 
dense regions are hotter than the mean density gas) after reionisa- 
tion. The temperature Tigm.o at reionisation is determined by the 
spectrum of the background radiation: it is higher the harder the 
spectrum. If the spectrum does not evolve afterwards, the tempera- 
ture of the IGM is expected to decrease almost adiabatically (e.g., 
(7 - 1) ~ 0.6). 

In our models the adopted spectrum at z rc i — 7 is consis- 
tent with the temperature and (7— 1) at £ = 3 — 5 derived from 
the line widths of the Ly man-q forest: Tigm.o ~ 1 — 3 x 10 4 
K, (7 - 1) ~ 0.1 - 0.4 iRicotti et alJhoOCtjschave et alJl200Ct 
ITheuns et ai1l2002h . We find that the IGM temperature decreases 
slowly with redshift (Tigm oc (1 + z) ' 5 ) and the equation of 
state has (7 — 1) = 0.3 ±0.1 (i.e., an intermediate value between 
the isothermal and adiabatic case). Loose constraints (due to the 



1 This is due to the shorter cooling time at higher densities. At Tigm ~ 
10 4 K, H2 cooling is dominant over the Lyman-a cooling. H2 cooling is 
not very important and the dependence of the temperature on the density is 
weak. 



large errors on the measured temperature) can be put on reionisa- 
tion models using the observed Tigm- If reionisation happens too 
early and the IGM is not reheated by additional sources that pro- 
duce the hardening of the background spectrum, the temperature 
at z ~ 3 — 5 would be too low when compared to observations. 
Or, if the spectrum is extremely hard at z la i ~ 7, the IGM tem- 
perature at z ~ 3 — 5 would be higher than observed. Compton 
heating is always much less efficient than photoionisation heating 
in these models. We have estimated that it is subdominant if the 
neutral fraction is xu 1 Si 10~ 7 . 

In Fig.|6|we show the He ionisation history for model M3. The 
dotted, solid and dashed lines refer to gas over-densities 8 — 0.1,1, 
and 10, respectively. Interestingly, in this model (and in model M4) 
the He II was almost fully reionised at high-z and fully reionised at 
z ~ 3. Helium double ionisation is almost complete at z ~ 15 
but afterwards Hell recombines very slowly as the X-ray emis- 
sivity decreases. Note that the recombination rate of He III is four 
times faster than the one of H II . The reason for this slow recom- 
bination rate of He III is the following. In this model the second 
reionisation of He II it is not produced by a second peak in the X- 
ray emissivity at z ~ 3 due to AGN activity. Instead it can be 
explained by the hardening of the spectrum of the background ra- 
diation in the UV bands, that it is produced by the combination 
of two effects: (i) the IGM becomes increasingly transparent to 
ionising radiation after Hi reionisation at z Ic i = 7 and (ii) the 
background soft X-rays are redshifted into the hard-UV bands and 
ionise more efficiently the He II . This effect is a general feature 
of all late and intermediate accretion models. This is illustrated in 
Fig-0 that shows the effect of the redshifted background on the 
evolution of the He I and Hell photoionisation rates, F. Neglect- 
ing recombinations and the Hubble expansion we have approxi- 
mately that (tn /n(He))dn(He)/dt ~ i^T, where tn is the Hub- 
ble time. The solid lines show £.HT(He Ii) for He II and the dashed 
lines tflT(Hel) for He I as a function of redshift. The thick lines 
refer to model M3 and the thin lines to model M4. 

Observations of the He II Lyman-a forest (e.g.. lReimers et all 
1 19971 ITheuns et al| 2002 ) and the line widths of the H I Lyman-a 
forest (Ricotti et al. 2000) at z ~ 3 suggest that Hell reionisation 
happens at redshift z — 3. This is usually attributed to the observed 
peak o f AGN activity at that redshift, but as noted in Mini ati et alJ 
(2004), heating from shocks arising from cosmic structure forma- 
tion also can make a significant contribution. The scenario pre- 
sented here is another available mechanism that can explain the 
observed He II reionisation at z ~ 3 and is consistent with the red- 
shift evolution of IGM equation of state. Note that in our models 
the temperature decreases slowly, monotonically, but does not show 
any substantial increase corresponding to the redshift of He II ion- 
isation at z ~ 3. The same smooth behaviour is exhibited by the 
parameter (7 — 1). 

The evolution of the SFR for models M3, M4 and M5 was 
shown in Fig. [2] Note that in model M5, in which BH accretion 
starts earlier, the thermal feedback strongly suppresses star forma- 
tion in the smaller mass galaxies. This effect reduces the number 
of seed BHs and therefore the emission of X-rays at high redshift. 
The filling factor of H II regions (shown in the top panels of Fig.[5} 
has an evolution similar to that of the SFR: after the initial growth, 
the filling factor remains constant at a few percent of the volume. It 
starts growing again when most massive galaxies form at z ~ 15. 
This effect is also observed in the cosmological simulations with 
radiative feedback discussed in paper lib. 

The optical depths to Thomson scattering produced by the X- 
ray preionisation models are within the 65% confidence limit of 
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Figure 6. Helium ionisation history for model M3. The dotted, solid and 
dashed lines refer to gas over-densities 6 = 0.1, 1, and 10, respectively. 
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Figure 7. Effect of the redshifted background on the evolution of the He I 
(dashed lines) and He II (solid lines) photoionisation rates, I\ The lines 
show (t[j /n(He))dn(He)/dt ~ iflT, where tu is the Hubble time, 
for model M3 (thick lines) and model M4 (thin lines). 



WMAP. The models are also in agreement with r c ~ 0.103±°'o4° 
found by Teemarketal. (2003) by combining WMAP and other 
CMB experiments with the Sloan Digital Sky Survey (see their Ta- 
ble 4). The optical depth, r c , and the visibility function g(z) = 
e~ T ° dr c / 'dr\ (where r\ is the conformal time) as a function of red- 
shift for models M3, M4, and M5 are shown in Fig. [H] The three 
models differ for the accretion history onto seed black holes. In 
model M5 most of mass in black holes is accreted at z > 25, in 



Figure 8. Thomson scattering optical depth, t c , and visibility function, 
g(z), as a function of redshift. Models M3 (solid lines), M4 (dotted lines) 
and M5 (dashed lines) are shown. The models differ for the accretion history 
of BH seeds (duty cycle) but have the same total integrated X-ray energy 
input. The maximum of r e is obtained for model M4, that has the bulk of 
accretion at intermediate redshift (z ~ 20) with respect to models M3 and 
M5. 



model M4 at z > 20 and in model M3 at z > 15. It is interest- 
ing that model M4 has the larger r e , even if all the models have the 
same integrated energy input from accreting BHs. This is due to the 
small, but non-negligible, recombination rate of partially ionised 
gas. If the IGM is partially ionised at high redshift and recombi- 
nations can be neglected, r e should increase the earlier the partial 
ionisation begins, and vice- versa if the partial ionisation starts later. 
But if recombinations are not negligible at high redshift, then there 
is an intermediate redshift that maximises r c . In our models this 
redshift is z ~ 15 (model M4). 

Another interesting effect worth noticing is the delay between 
the formation of the first sources and the build-up of a substantial 
background. Roughly the time scale for building the background 
is given by a fraction of the Hubble time. For instance, if the first 
sources formed at z = 40(30), the X-ray background will build up 
at z — 25(20). Also, as already noted, the age of the universe at 
z — 30 is about one Eddington time t^d = 10 8 yr. If seed BHs 
accrete at the Eddington limit their mass would grow exponentially 
one e-folding from z = 30 to z = 20 and 4 e-foldings to z = 
12. Therefore, the energy available to produce X-rays in the very 
early universe is limited to 5-10 times the mass of the seed BHs. 
This means that we need an increasingly large mass in seed BHs 
to preionise the IGM starting at higher redshift (cf., Fig.|5J- This is 
only possible if the IMF of the first stars is top-heavy and has been 
discussed in paper I. 

To illustrate the effect of cosmological parameters on the r 
produced by X-ray preionisation, in Fig.|5|we show r c and g(z) for 
model M3b that has a spectral index of initial density perturbations 
n s — 1.04 and, for comparison, model M3 that has n s — 1 nor- 
malised to the same amplitude o f the power spectr um at k = 0.05 
h Mpc -1 measured by WMAP l iVerde et alJl2003l) . The other pa- 
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Figure 9. As in Fig. |8} for model M3 with spectral index of primordial 
perturbations ra s = 1 (solid lines) and the same model but using n B = 1.04 
(model M3b, dot-dashed lines). 



Table 3. Results of the semianalytic models 



Model 


Te 


max 
Z 9 


y-par 


^BH 


/XRB 


r(HeIl)t H 








xlO 7 


(%) 


(%) 


(%) 


M3 


0.146 


12.0 


3.64 


3 


5(12) 


35 


M3b 


0.166 


12.5 


4.26 


3 


5(12) 


35 


M4 


0.170 


17.0 


3.95 


9 


7(17) 


25 


M4b 


0.170 


17.0 


3.95 


9 


7(17) 


25 


M5 


0.140 


21.5 


2.96 


42 


0.3(1) 


1.5 



Meaning of the values in each column: r c is the Thomson scattering op- 
tical depth; z™ ax is the high-redshift maximum of the visibility function 
g(z); y is the Compton distortion parameter; wbh/ w * is the maximum ra- 
tio of BH to stellar cosmic density (to translate i^bh to PBH multiply by 
5.5 X 10 9 Mq Mpc — 3); /xrb is the fraction of the X-ray background 
at 50-100 keV (and 2-10 keV, in parenthesis) due to early black holes; 
r(HeIl)tH ~ (tH /nHeu)dnii su /dt is the fractional rate of Hell pho- 
toionisation per Hubble time per helium atom due to the redshifted X-ray 
background. 

rameters in model M3b and M3 are the same. In model M3b seed 
BHs form earlier and the preionisation starts at higher redshift than 
in model M3. In this model we get r c = 0.16, compared to r c 
= 0.14 of model M3. This result is qualitatively different if reion- 
isation is instead produced by UV sources. In this second case, in 
paper I we have shown that r c is not substantially increased if the 
power spectrum of initial perturbations has more power on small 
scales because star formation in small mass galaxies is regulated 
and limited by feedback effects. 

In conclusion, adopting the model constraints listed at the 
beginning of this section, we obtain values of r c consistent with 
WMAP. The models are consistent with and do not violate any ob- 
servation at z < 6 that are evident. In the following paragraphs we 
discuss a few obvious observational constraints. But we will focus 
and return to this discussion in paper lib. 

Reheating of the IGM at redshifts z < 10 5 , produces a Comp- 
tonised spectrum of the CMB blackbody. The Compton distortion 



Figure 10. Background radiation at z = 12, 6, 3 and 2 = for model 
M4. The thick line shows the X-ray background radiation at z = 
IComastri et alll995l) . The redshifted X-ray background produced by early 
epoch of black hole accretion constitutes about 5 — 10% of the hard X-ray 
background in the 2-50 keV bands. 



is usually parametrised by the value of the y: 

_ f' k B (T e — ToMB), , -7^, n -6 f T " T e 
V = J ^ dr e ^ 1.7X10 ^ —dr e , 

where T e is the electron temperature, Tcmb is the CMB photon 
temperature and r e is the optical depth to electron Compton scatter- 
ing. An upper limit on the y-parameter y 1.5 x 10~ 5 at 95% CL 
has been determined using COBE satellite data jFixsen et alll99(j) . 
The maximum value of y that we find in our models, i/~4x 10~ 7 , 
is not constrained by the current data on the spectral distortion. A 
future experiment with an upper limit on the y-parameter improved 
by two orders of magnitude would be able to detect a possible X- 
ray preionisation and reheating in the early universe. 

InFig.Eol we show the background radiation at z = 12, 6, 3 
and z = for model M4. The thick line shows the X-ray back- 
ground radiation at z — iComastri et alj|l995T) . We see that the 
redshifted X-ray background produced by early epoch of black hole 
accretion constitutes about 5 — 10% of the hard X-ray background 
in the 2-50 keV bands. 

The observed (X,7) ray background is produced almost com- 
pletely by known sources (mainly Seyfert galaxies at z ~ 1 — 2). 
Howe ver there is some roo m for additional sources and recently 
iDe Luca & Molendil j2004l) have estimated that that about 20% of 
the observed background may be produced by a new population 
of faint X-ray sources, currently undetected within the sensitivity 
limits of the deepest X-ray surveys. The intensity of the high red- 
shift X-ray background decreases from its value at redshift z as 
(l+z)~ 3+< *, where (f> is the logarithmic slope of the spectral energy 
distribution. By imposing that the redshifted background produced 
at the redshift of the peak of the BH accretion does not exceed 20% 
of the observed value at z — we find the following upper limit 
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for the BH accretion rate: 

p-r Z (lO- 2 M Mpc- V" 1 ) (y ) ~ 3 ' 7 • (ID 

Note that we this expression has been derived for the accretion his- 
tory of the models presented in the present paper. It implicitly as- 
sumes that p ac = A rho^/ At w ' m At ~ tij(z), where tii(z) 
is the Hubble time at the redshift of maximum accretion (i.e., the 
accretion rate rapidly fades after the maximum). 

In table|3|we summarise the values of some relevant quantities 
computed from the semianalytic models listed in table |2| See the 
table footnote for the explanation of the meaning of the quantities 
shown in each column of the table. 



5 DISCUSSION AND SUMMARY 

In a previous paper (paper I), we showed that even with a top heavy 
mass function Population III stars could not reionise the universe to 
a significant degree if they end their lives as pair instability super- 
novae. The reason is that such events contaminate the high density 
regions to such an extent that further formation of extremely low 
metallicity stars becomes impossible before significant ionisation 
has occurred. But, if these same stars implode to massive black 
holes, then accretion onto them at near Eddington rates will pro- 
duce a pervasive X-ray background which can quite efficiently par- 
tially preionise the bulk of the cosmic baryons before z — 15 and 
give an optical depth to Thompson scattering (cf., table [5} consis- 
tent with WMAP results. The primary reason for the increased effi- 
ciency is the long mean free path of X-rays which both increases the 
fraction of photons that escape the parent halo and ionises preferen- 
tially in the low density regions where recombination rates are low. 
In addition X-rays only partially ionise the IGM to 10-20%, there- 
fore the recombination rate is reduced with respect to the reion- 
isation scenario produced by UV sources. Thus we find that the 
average number of times a baryon is photoionised in this picture is 
about unity for r e ~ 0. 15 while it is about ten in the case when stel- 
lar UV is used (paper I). Moreover, since the X-rays preionise the 
low density regions of the IGM at early times, the required number 
of ionising photons from Population II stars to have reionisation at 
redshift z — 6 is much reduced. 

But we might ask, could the X-ray s arise as e asily from X- 
ray binaries or could SN remnants (cf., Oh 2001)? These other 
sources are certainly helpful but they will generate somewhat 
less energy in X-rays than accretion can do. The ratio of X-ray 
energy to massive star formation expressed as an efficiency is 
£SN = -B 5 ilO 51 erg/143M c 2 ~4x 10~ 6 for these sources 
(assuming Population II with Salpeter IMF and neglecting Comp- 
ton cooling of the X-ray emitting hot gas) but we see from e ac = 
0.15/3(MsMBH/A4bui gc ) that the comparable number from accre- 
tion is 2.25 x 10" 5 . 

What are the observable consequences of the scenario we are 
proposing (in addition to generating a large r c )? 

(i) The power spectrum of the cross correlation of temperature 
and polarisation of the CMB is sensitive to the visibility function 
g(z) = e _Tc dr c / 'drj (cf., Fig. l8l9i . The Planck satellite should be 
able to distinguish between the visibility function produced by an 
early X-ray partial ionisation (peaking at large z) or the one ex- 
pected for reionisation by stellar sources (peaking at lower z). 

(ii) Since the assumed spectrum (Fig. [3) has a positive slope 
(d\aL u /dhxv) > 3 for much of the relevant range of frequen- 
cies there are ever fresh ionising photons available and the proba- 



bility per unit time that a neutral hydrogen atom is ionised does not 
decrease with time. This produce softer photons that, even if the 
ionisation efficiency of X-rays decreases when the ionisation frac- 
tion exceeds 10-20%, enhance the ionisation fraction in the voids 
to values larger than 20 % because of the contribution from the 
redshifted photons of the X-ray background. In most models, the 
photons from the redshifted X-ray background can fully reionise 
He II at redshift z ~ 3 without any additional contribution from 
quasars at lower redshifts. Due to the heating rate by He II ionisa- 
tion the temperature of the mean density intergalactic medium re- 
mains close to 10 4 K down to redshift z ~ 1 and the IGM equation 
of state has index (7 — 1) ~ 0.3. 

(iii) The maximum accretion rate onto seed black holes as a 
function of redshift is limited by the observed magnitude of the 
7 ray background at z = 0. This limits the ability of late preioni- 
sation models to produce large values of r e and contribute substan- 
tially to the mass growth via accretion of seed black holes. In mod- 
els that produce optical depth r c ~ 0.17 consistent with WMAP, 
the redshifted X-ray background produced by the early epoch of 
black hole accretion constitutes about 5 — 10% of the hard X-ray 
background in the 2-50 keV bands. 

(iv) The model predicts that dwarf spheroidal galaxies, if they 
are preserved fossils of the first galaxies, would host a mass in 
black holes that is 5-40% of their stellar mass (cf., table E). The 
larger values are produced in the early preionisation scenarios and 
it might be possible to rule them out on the base of observations 
from dynamical considerations. 

(v) The IGM is reheated to T = 10 4 K before the ionisation 
fraction exceeds 10% producing typical values of the Compton dis- 
tortion parameter t/~4x 10~ 7 . The observational upper limit on 
this parameter is y < 1.5 x 10~ J , too large to constrain our models. 
Star formation in the smaller mass haloes is reduced by the increase 
in the IGM Jeans mass following this early reheating. 

(vi) The redshifted 21cm signal in emission and absorption 
against the CMB can be used to discriminate between preionisation 
scenarios and reionisation from UV sources. This will be discussed 
in paper lib, where we will also show that preionisation by X-rays 
produces CMB secondary anisotropics on small angular scales that 
are easy to recognise when compared to models of inhomogeneous 
early reionisation by stellar sources with the same r c . 

What are the requirements of the X-ray preionisation models 
in terms of seed black hole production from the first stars and ac- 
cretion? 

Since the time scale for accretion, £e<i = 10 s yr, is compara- 
ble to the Hubble time at 2 = 30 the maximum X-ray emissivity 
at high-redshift is proportional to the total mass of the seed BHs. 
For this reason, in order to get a larger r c in early preionisation sce- 
narios, a top-heavy IMF is favoured over a Salpeter IMF. But this 
is true only if super-massive stars mostly end their lives collapsing 
directly into BHs without exploding as SNe or as pair-instability 
SNe. 

We have shown models where the accretion onto seed black 
holes is near the Eddington rate at early redshift. This might per- 
haps be difficult to achieve. But it is worth noticing that it is not 
strictly necessary that the accretion is at the Eddington rate initially. 
Although, if the accretion rate is less efficient than that, a larger 
stellar mass fraction producing seed black holes, /bh, is needed to 
accommodate for this (see equation {SJ). 

In this paper we have explored a range of models with dif- 
ferent accretion histories onto seed black holes. But we note that 
the models are not all equally plausible scenarios. The very early 
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preionisation scenario (M5 ) is not the most f avourable for several 
reasons: (i) as noted by Madau et alj 120041) . accretion in small- 
mass haloes can be difficult because of the substantial gas photoe- 
vaporation in these galaxies (ii) if seed BHs do not accrete at near 
the Eddington rate the mass fraction, /bh, of seed BHs is extreme 
also assuming top-heavy Population III stars (iii) the number of re- 
combinations per ionising photon in this model is larger than in 
the intermediate preionisation model and therefore r e is reduced. 
The late preionisation model (M3) is also not favourable because 
it produces a low r c if we constrain the maximum black hole ac- 
cretion rate using equation i 1 1 i to be consistent with the observed 
7-ray background at z = 0. We therefore conclude that preionisa- 
tion starting at z ~ 15 — 20 (e.g., model M4) is the most plausible 
scenario to reproduce r c measured by WMAP. In this case ~ 10% 
of the X-ray background in the 2-50 keV band is produced by the 
early generation of mini-quasars and roughtly half of the currently 
estimated black hole mass density was formed at early times (even 
if only a fraction of them may have merged into the SMBHs in the 
galactic centres). In paper lib we study this scenario using hydro- 
dynamic cosmological simulations with radiative transfer and we 
discuss observations that could probe this model. 

We also note in conclusion that this picture while consistent 
with the central values of WMAP results could not easily produce 
T e at the upper end of the presently allowed values. If r c Si 0.2, 
X-ray preionisation is not likely the dominant responsible mech- 
anism and other scenarios need to be considered (e.g., decaying 
particles, non Gaussian perturbations, spectral index n s > 1, pri- 
mordial BHs, unconventional recombination). 
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APPENDIX A: A SEMIANALYTIC MODEL FOR 
REIONISATION 

We implemented a semianalytic model to study reionisation, chem- 
ical evoluti on and re-heating of th e IGM. The code is based on the 
method of IChiu & OstrikeJ i20obT) for the calculation of the fill- 
ing factor of ionised regions and for the star formation recipe. The 
main difference is the inclusion of radiative transfer for the back- 
ground radiation that allows us to follow the thermal and chemi- 
cal evolution of the IGM outside the H II regions surrounding each 
UV source. We also consider the thermal feedback produced by 
the X-ray reheating of the IGM (see S IA4> . The mass function of 
DM haloes and their formation/merger rates are calculated using 
the extended Press-Schechter formalism. The star formation rate 
is assumed to be proportional to the formation rate of haloes (see 
S IAli . We consider the IGM as a two-phase medium: one phase is 
the ionised gas inside the H II regions and the other is the neutral 
or partially ionised gas outside. We solve the radiative transfer for 



the volume averaged specific intensity and we derive the specific 
intensity inside and outside the H II regions from their volume fill- 
ing factor and by separating the contribution from local UV sources 
and the background radiation from distant sources (see §|A3j. We 
solve the time-dependent chemical network for eight ions and the 
thermal evolution as a function of the gas density outside the H II 
regions. Given the clumping factor of fully ionised gas around the 
UV sources we evolve the filling factor, temperature and chemistry 
inside the H II regions (see S IA2I . The following heating and cool- 
ing processes are included: collisional- and photo-ionisation heat- 
ing, ionisation by secondary electrons (see § IA2.lt . H, He, H2 cool- 
ing, Compton and adiabatic cosmological cooling. The rates are the 
same as the ones used in tRico tti et all200l1) . For some test cases 
the results of the semianalytic code are in good agreement with the 
results of the cosmological simulations with radiative transfer pre- 
sented in this paper and in lRicotti et alJl2bo2bh . 



Al Star Formation 

We assume that stars and quasars form in virialised DM haloes. 
We use the extended Press-Schechter formalism to calculate the 
formation rate of bound DM haloes of a given mass Mdm at the 
time t. In general, given a probability distribution function (PDF) 
of density perturbations, P(y), the DM comoving mass fraction of 
collapsed objects with mass between M dm and Mdm + dM dln is 
given by Q, ps dlnM dnl = P(y)du, where 



d lnMdm^ps 



dvP(v) 



The comoving number density of collapsed objects N ps = 
^psPo/Afdm, where po is the mean DM density at z — 0, is given 
by 



N pB dM dm 
V 



A (Aid 



2p P(v)vA(M dm _)dlnM da 

a{M dm )D(t) ' 

dlnv dlno-(Mdm) 
dlnM dm ~ 



din M d 



(Al) 
(A2) 

(A3) 



where a(M dm ) is the variance of the fluctuation on a scale M dln 
linearly extrapolated to z = 0, 8 C ~ 1.69 is the linear overdensity 
for collapse of a top-hat perturbation and D(t) is the linear growth 
factor. Here we use a Gaussian PD F P = (2ty)~ 1 ^ 2 exp (— u 2 /2), 
but the method has been applied by IChiu & Ostriker 2000) to non- 
Gaussian PDFs. With a similar procedure as in lRlcottl2002l) . we 
normalise the star formation efficiency by assuming that at z — 
a fraction of baryons w» = 14% has been converted into stars. 
The global SFR between redshift < z < 5 agrees within the 
errors with the observed global SFR (e.g., Lanzetta et al. 2002J). As 
shown in §0 a similar method of normalisation is applied to the 
BH accretion efficiency by assuming that at z — a fraction of 
baryons ubh = 10 -4 is in black holes. 

The formation rate of virialised objects, Nf orm , can be de- 
rived using the Press-Schechter formalism and calculating the rate 
of "destruction" of bound objects that are incorporated in larger 
haloes: 



= N ps + N d 



(A4) 



As in IChiu & OstrikeJ l2000l) . we use a destruction probability 
4>(t) = N dcEt /N ps = D/D derived assuming that the destruc- 
tion probability is scale-invariant (see Sasa kill994l) . It follows that 
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the probability that an object formed at time t / exists at time t is 

rt! D{t) 



pBUTv(t\t{) = J ' dt<t>(t) = 



D{tiY 



(A5) 



where D(t) is the linear growth factor. At redshifts z > 1, 
Psurv(i|tf) = a(t)/a(tf) for the concordance ACDM cosmology. 
From equation <A4t . recalling that D/D « f2°' 6 (d/a), it follows 

N !om = -N ps ^4>(t) = N ps u 2 n ' 6 (- 
P dv \a y 

where d\nP/du — —v for a Gaussian PDF (P oc exp (—v 2 /2)). 
The comoving number density of haloes at time t that formed at 
time t f is given by 

N(Mdm,tf,t)dMdmdtt = N{ OI m(Mdm,tf)p B ui V (t\t{)dMdmdt{. 

The emission rate per unit volume from the sources, often called the 
source function S(t), when expressed in dimensionless units S = 
S{i) / (nnhvo), where hvo = 13.6 eV and uh is the hydrogen 
mean number density, is given by 



{■ OC {' t 

S(t) = no 1 / dMdm / dtt 
Jo Jo 



(A6) 



x*(M dm ,t f )iV((M dm ,t / ,t)£(M dm ,tf,t). 

The function $(Afd m ,tf), defined in S IA4I is a step function that 
determines the minimum mass of forming galaxies as a function of 
the their time of formation, tf. The dimensionless luminosity from 
each source is C v oc g „ €\jy M*c 2 x exp (— x/tdyn), where g v is the 
SED of the source, x = (t — tf)/tdyn and tdyn = (37r/32Gp v ir) 
is the free-fall dynamical time. 



A2 Cooling and Chemistry 

We solve the time-dependent equations for the photo-chemical for- 
mation/destruction of eight chemical species (H, H + , H~, H2, Hj , 
He, He + , He ++ ), inclu ding the 37 main process es relevant to deter- 
mine their abundances I Shapiro & Kana I l987h . We use ionisation 
cross sections from Hui & Gnedhj 119971) and photo-dissociation 
cross sections from A bel et alJll997l) . We solve the energy conser- 
vation equation 



dE e . 



dt 



= r-A, 



(A7) 



where E gas = (3kT/2)n H (l + x(He) + x e ), with n H the hy- 
drogen number density, and x(He) and x e the helium and elec- 
tron fractions, respectively. Note that uh is a function of time be- 
cause of Hubble expansion. Also, x e is time-dependent; neglect- 
ing this effect results in a temperature that is overestimated by 
about a factor of two. The cooli ng function, A, include s H and 
He line and continuum cooling IShapiro & Kandll987h . H2 ro- 
tational and vibrational cooling excited by collisions with H and 
H 2 iMartin et alJll996l : iGalli & Pall al Fl998> and adiabatic cosmic 
expansion cooling. The heating term, V, includes Compton heat- 
ing/cooling and photoionisation/dissociation heating. We solve the 
system of ODEs for the abundances and energy equations, using a 
4 th -order Ru nge-Kutta solver. We switch to a semi-implicit solver 
(Gnedin & Gnedinl ll998l) when it is more efficient (i.e., when the 
abundances in the grid are close to their equilibrium values). The 
spectral range of the radiation is between 0.7 eV and 10 keV. The 
primordial helium mass fraction is Yp = p(He) /pt, = 0.24, where 
Pb is the baryon density, so that x(He) — Vp/4(1 — Yp) = 
0.0789. The initial values at 2 = Zi for the temperature and species 



abundances in the IGM are: T = 10 K, xh 2 = 2 x 10~ 6 , and 
x e ~ x H + = W 5 /(hn b nl /2 ) = 6.73 x 10" 4 . The initial abun- 
dance of the other ions is set to zero. 

Outside the Stromgren spheres surrounding the UV sources 
the thermal and chemical histories are calculated as a function of 
the gas density, using the specific intensity (J„)back of the back- 
ground radiation defined in § IA3.1I In this paper we show the abun- 
dances and the temperature of the IGM for three values of the 
baryon overdensity S — p/po = 0.1, 1, 10. Inside the Hll re- 
gions we calculate the abundances of only H and He ions, since 
the abundances of molecular hydrogen and its ions are negligible, 
being photodissociated by H I ionising radiation. The effective co- 
efficients of collisional ionisation and recombinations are larger by 
a factor Chu, to take into account the gas clumping inside the H II 
regions (see S IA3I . The photoionisation and photoheating rates are 
calculated as shown in § IA2.1l using the specific intensity Jhh of 
the radiation field inside Hll regions. The clumping factor Chu is 
tabulated as a function of redshift by fitting the clumping in the 
cosmological simulations with radiative transfer presented in pa- 
per lib. 

A2. 1 Secondary Ionisation and Heating from X-rays 

Photoionisation of H I , He I , and He II by X-rays and EUV photons 
produces energetic photoelectrons that can excite and ionise atoms 
before their energy is thermalised. Th is effect can b e important be- 
fore reionisation lohll200lllVenkatesan et all200ll) . when the gas 
is almost neutral and the spectrum of the background radiation is 
hard due to the large optical depth of the IGM to UV photons. 

Collisional ionisation and excitation of He II by primary elec- 
trons are neglected since, in a predominantly neutral medium, 
they are unimportant. The primary ionisation rate for the species 
i = Hi, He I, He II is, 



c 



4ty 



J, 



exp(—T„)aldv, 



(A8) 



where r„ is the continuum optical depth, and al is the photoioni- 
sation cross section of the species i. Secondary ionisation enhances 
the photoionisation rates as follows: 

C H, = C H '+ E C(^\Elxe)) (A9) 

i=Hl,HeI,HelI 

C He. =c Ha, + £ i^' (E^, X e )) , (A10) 

i— Hi,Hei,He n 

where {$ Hl (^,a; e )) and {$ HcI (^o, x e )) express the average 
number of secondary ionisation per primary electron of energy 
Eq = h v v — I 1 weighted by the function (J v /h p u) exp(—r v )al. 
Here P — h v Vi is the ionisation potential for the species i. 

The photoionisation heating rates for the species i — H I, He I 
are given by, 



r l = 4tt 



/ -^exp(-T)alE h {Eh,Xe)dv. (All) 



Analytic fits to th e functions $ Hl , $ Hel and Eh, based on the Monte 



Carlo resultsj)f JSlrull&_TOn_Steenberj 1 19851) are given in Ap- 
pendix B in lRicotti et alj <2002ah ! 



A3 Radiative transfer 

The evolution of the specific intensity J„ [erg cm -3 s _1 Hz" 1 sr" 1 
] of ionising or dissociating radiation in the expanding universe, 
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with no scattering, is given by the following equation: 



dJ v d ( .i \ ( dJ v 



- hi/ J v -\- Sv 



(A12) 

Here, x 1 are the comoving coordinates, H is the Hubble constant, 
k v is the absorption coefficient, S v is the source function, and 
x 1 — en} I a, where n l is the unit vector in the direction of photon 
propagation and a — (l + ^) _1 is the scale factor. The volume- 
averaged mean specific intensity is 



J„(t) = (J v (t,x,n)) v , 



(A13) 



where the averaging operator acting on a function f(x, n) of posi- 
tion and direction is defined as: 

(f(x,n)) v = lim / d 3 x [ dnf(x,ft). (A14) 

The mean intensity J v (i) satisfies the following equation: 

^-H^^-iJ^j = -k v J v +S v , (A15) 

where, by definition, S v = (Sv)y, and k v = (k v J v ) v / J v - In 
general, k v is not a space average of k v , since it is weighted by the 
local value of the specific intensity J v . But when the mean free path 
of radiation at frequency v is much larger than the characteristic 
scale of the problem (in our case the size of the computational box), 
we have k v = (k u )v- 

If we rewrite equation IA15i in terms of the dimensionless 
comoving photon number density n„ = (aLbox) 3 47rJ„/fo p where 
Lbox is a spatial scale (in this paper we take it to be the comoving 
size of the computational box), and h p is the Planck constant, using 
substitutions dt — (a 2 / Ho)dr and £ = ln(^), equation <A15> can 
be reduced to the following dimensionless equation: 



dne 1 da dnc 
or a dr at, 



(A16) 



where — a 2 k v /Ho and S% = a 5 L^ ^4irS v / h p Ho. Using a 
comoving logarithmic frequency variable, 



f=£ + ln(a), 



equation <A16> can be reduced to 



— — a^n^ + S( 



which has the formal solution, 



n^(r + At) — n^(r) exp 



+ 



dt'S^(t') exp 



dt'a0) 
dr'a^ (t') 



(A17) 



(A 18) 



(A19) 



We calculate equation <A19> at each time step, At, of the simula- 
tion. The two integrals inside the square brackets on the right side 
of equation <A19> can be solved analytically. We solve the third 
integral numerically. 

A3.1 Background radiation and radiation inside H II regions 

Two terms contribute to the volume averaged specific intensity in 
equation <A19> . The first term is produced by the background ra- 
diation from redshifted distant sources and the second term is pro- 
duced by the local sources: 



t Tsour . jback 



(A20) 



If the emissivity of the sources is zero (<S = 0), it follows that 
jsour _ q Q m y pjj t ons with a mean free path larger than the 
mean distance between the sources contribute to the background 
radiation, JJ ack . 

If the volume filling factor of the Stromgren spheres surround- 
ing the UV source is /mi, we can rewrite the volume averaged spe- 
cific intensity as the sum of the mean specific intensity inside the 
volume occupied by the Hll regions and the volume outside the 
H II regions: 



JtH II r I rback/-] 
v = Jv /Hi] + J v (1 



/HI.). 



From equation JA201 and equation JA2U we get, 

J "sour 



TH I 

J J J 



/hii 



. yback 



(A21) 



(A22) 



In order to derive the specific intensity, J„ ", inside the Stromgren 
spheres we need to estimate their volume filling factor /hh. We 
calculate the evolution of the fillin g factor of H II regions f ollowing 
closely the method introduced by Chiu & OstrikerlT200fj) . For the 
sake of completeness we show the equations that we solve to derive 
/hh in the next section. 



A3. 2 Filling factor of H II regions 

In order to calculate /hh it is more convenient to introduce the 
porosity parameter, Q(t), of the Hll regions defined by the rela- 
tionship /hh — 1 — exp [~Q(t)]. The porosity Q(t) offers the 
advantage that when /mi <S 1 we have Q(t) « /mi and when 
/mi — 1 we have Q(t) ^ 1. The porosity of the Hll regions is 
given by 



Q(t) 



dM dm / dt f 

u JO 



(A23) 



x^{M Am ,tt)N{{M dm ,t } ,t)VMn{M Am ,tt,t), 



where Vhh = (47r/3)r mi is the comoving volume filled by the 
Stromgren spheres, where run is their comoving radius. Compar- 
ing equation JA6> to equation <A23> we find the following relation- 
ship between the ionised volume per unit luminosity and the mean 
comoving luminosity density S: 



n Vmi _ Q(t) 

c(t) - s{t) ■ 



(A24) 



The contribution of each source to the mean specific energy density, 
E = A-nJ v /c is given by, 



c 



dM dm / dtf^(M dm ,tf)N(M dm> t f ,t) 



dVL 



dM dn 



47r(ar) 2 

dt f JV(M dm , U ,t)Lr Hu a(t)- 2 . 

(A25) 



Using equation <A24t we find that, in dimensionless units £ — 
E/(n,H hvo), the mean UV energy density is 

1/3 



a (3Q(t)\ 
c \ 4tt5 



-4/3 



dM dm / dt{ 

Jo (A26) 

4/3, 



x*(M dm ,t f )Ar(M dm ,^,t)£ 4 ^(M dm ,t f ,t). 

The time evolution of the mean UV energy density is calculated 
solving the energy conservation equation, 



d£ 
~dt 



7- + XHiCo- e n H 
a 



dfu 
dt 



(A27) 



© 2002 RAS, MNRAS OOO mHTl 



16 M. Ricotti and J. P. 'striker 



We have assumed that the spectral energy distribution of the ion- 
ising radiation is g v oc j/~ 7 . Here <r e is the effective UV energy 
loss cross section, defined as o e — (r p h + ri on )/(47r J). The time 
derivative of the filling factor is fn = (1 — fn)Q(t). The porosity 
parameter Q(t) and the volume filling factor of the H II regions are 
obtained solving equation <A27t and equation <A26> . 



A4 Radiative Feedback 

Feedback mechanisms on star formation are the most difficult part 
to implement using the semianalytic approach. The recipes used 
should be based on results of numerical simulations or scaling rela- 
tions based on observations (such as the dependence of the mass to 
li ght ratio on th e luminosity of galaxies). Analysing the simulations 
in lRicotti et alJ (2002b), it appears that star formation is regulated 
by the interplay of several processes: 

(i) Internal feedback by UV radiation that produces photo- 
evaporative winds 

(ii) Internal feedback by SN explosions 

(iii) External (but local) feedback from UV radiation that regu- 
lates H2 formation/destruction and cooling 

(iv) Thermal evolution of the IGM: the IGM Jeans mass sets the 
minimum mass of galaxies forming at a given redshift. 

In the semianalytic code, due to impossibility of modelling 
local external feedback processes, we only implement process (iv) 
which is global. We neglect the effects of internal feedbacks (i) and 
(ii). 

Galaxies with mass smaller than the Jeans mass in the IGM 
cannot virialise because the pressure of the IGM prevents the gas 
from falling into their halo potential well. Therefore the formation 
rate of galaxies with masses Mdm < 10 8 Mq is suppressed if the 
IGM is reionised or heated by X-ray or hard-UV background to 
T ~ 10, 000 K. This process works together with H2 destruc- 
tion to suppress the formation of small-halo galaxies. If we take 
into account the finite time required for pressure to influence the 
gas distribution in the expanding universe, then the filtering mass, 
Mf,igm, wmcn depends on the full thermal history of the IGM, 
provides a better fit to the simulation results than the Jeans mass, 
which instead depends on instantaneous values of the sound speed 
(Gnedin 2000b). The filtering mass of the IGM is simply related to 
the Jeans mass by the relationship, 



M 2/3 



da'Mj /3 {a) 



l-l a - 

a 



1/2 



The Jeans mass is given by 



4-7T , ■) 

Mj = —p\ 3 j, 



(A28) 



(A29) 



where Xj = 2'KCstg is the Jeans length, iff = (iTvCp)" 1 ^ 2 is the 
free-fall dynamical time and c s is the IGM sound speed. 

We implement the feedback imposing a minimum galaxy mass 
M m in — Mp,iGM(tf) and computing the new SFR and the IGM 
temperature iteratively. Using a geometric mean for the temperature 
we make sure to achieve convergence. In practice we use a step- 
function kernel, 



<E>(M dm ,i f ) = 



0, if M dm < AfeyGM(if 



1, if Mdm ^ M F ,IGM(t{), 

any time we integrate over the mass function of DM haloes. 



(A30) 
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